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Abstract 

We discuss the equation of motion of the driven pendulum and generalize it to arbitrary driving 
angle. The pendulum will oscillate about a stable angle other than straight down if the drive am- 
plitude and frequency are large enough for a given drive angle. The emphasis is on the parameters 
associated with a simply made demonstration apparatus. 
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I. INTRODUCTION 



The general theory of the driven inverted pendulum with the drive angle restricted to 180° 
is given in Ref. l| and has been discussed in Refs. OQQ; the latter reference includes a good 
introduction to the physical system and further citations. In the following, the behavior of 
the harmonically driven pendulum will be described for any drive angle from the vertical. 

An inverted pendulum demonstration^ that is designed to be clamped to a table top 
becomes more interesting when hand held. We will examine the stability of this system 
as a function of the angular frequency, the drive amplitude, and the drive angle from the 
vertical. This type of demonstration is best used in a junior level classical mechanics course 
when introducing Lagrange's equations to find the equations of motion. Examples of driv- 
ing the support of a simple pendulum harmonically can be introduced, but only to set up 
the equation of motion. 6,7 The same apparatus can be used in an advanced graduate clas- 
sical mechanics course where the harmonically driven pendulum is used as an example of 
Lagrange's equations and as a source of problems.— 

While holding the saw in my hand,— with the power still on, I lowered the saw and 
observed the peculiar behavior of the pendulum as it sometimes found new stable angles 
of oscillation as the drive angle changed. Changing the driving angle in the plane of the 
pendulum's motion introduces the new and interesting problem that I address in this paper. 



II. EQUATION OF MOTION OF THE PENDULUM DRIVEN AT ANY ANGLE 

The dynamics of a driven inverted pendulum near 9 = 180° are described in Refs. QJ^QQ- 
We now generalize the problem to arbitrary driving angles. The geometry of the problem 
is illustrated in Fig. ^ A thin uniform rod of length L is driven at one end (pivot) with 
amplitude A, angular frequency ou, at an angle of 9d from the downward vertical. The angle 
of the thin rod is measured by the generalized coordinate 9 which also is measured from 
the vertical downward. (We also could have used a simple pendulum consisting of a mass 
m at the end of a light rod of length £, or a general pendulum with moment of inertia / 
and center of mass at a distance z from the pivot. The rod pendulum geometry chosen here 
corresponds to a simple lecture apparatus.^ A different pendulum geometry can be easily 
substituted into the equations of motion below.) 
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The angle of the rod 9 could be measured from the driving direction, 9 d . However, for 
only two values of 9 d (0 and tt) will the equilibrium angle of the pendulum be da- For any 
value of 9 d , there is a stable point for small uj and A at 9 = 0. So we choose to measure the 
pendulum angle 9 from the downward vertical. 

The kinetic energy and potential energy of the rod can be written in terms of the general- 
ized coordinate 9. The energies are separated into the motion of the center of mass plus the 
rotation about the center of mass. We first express the Cartesian coordinates and velocities 
of the center of mass in terms of 9. 

L , 
£ C m = +— sin 9 + A sin 9 d cosut (la) 

y cm = — — cos 9 — A cos 9 d cos uot (lb) 

^cm = 9 cos 9 — Aw sin 9 d sin ut (lc) 

Vcm = +2" 9 sin 9 + Alu cos 9 d sin cut. (Id) 

We can express the kinetic energy and potential energy in terms of the generalized coor- 
dinate 9. 

T = \m(x 2 cm +yl m )+ l -I c J 2 . (2) 
If we use the Cartesian coordinates and velocities from Eq. (0), we obtain 

X TTtL^ 

T = — [ 9 2 + mLAu9(sm 9 cos 9 d — cos 9 sin 9 d ) sin out + mA 2 uj 2 sin 2 cut] . (3) 

2 3 

The first term on the right-side of Eq. Q corresponds to the rotation about the moving 
support at the end of the pendulum with a moment of inertia of mL 2 /3 about that end. 
The potential energy depends only on y cm , and from Eq. ((TJ) we have 

V = mgy cm = 1— cos 9 — mgA cos 9 d cos tut. (4) 

The Lagrangian is £ = T — V, and using Eqs. Q and (jlj), we have 

£ = - [ 9 2 + mLAu9 sin (9 - 9 d ) sin ut + mA 2 u 2 sin 2 ut] 

2 3 

+ mg[— cos 9 + A cos9 d cos cot) . (5) 

The equation of motion is found from Lagrange's equation for the single generalized coordi- 
nate 9: 

d ,d£, dC 



If we evaluate the derivatives of the Lagrangian from Eq. (jSJ), we obtain the equation of 
motion for a pendulum driven at any angle: 

9 + sin(# - 6 d ) cos tot + ^ sin 9 = 0. (7) 

Note that if 9^ = 180°, we obtain the usual driven inverted pendulum . 1 i 2 i 3 i 4 This special case 
is revisited briefly in Sec. IIV Bl and in the Appendix. 

III. EFFECTIVE POTENTIAL 

We now introduce an effective potential to help us understand the physical origin of the 
stability of the inverted pendulum and to parameterize the condition for stability at any 
driving angle. 

Landau and Lifshitz 1 - separated the motion of the horizontally or vertically driven pen- 
dulum into two parts: a fast component at the drive angular frequency to, and a slow 
component <j)(t) that describes the slower overall swinging of the driven pendulum: 

9{t) = <j>(t) + £(t) (8) 

The angle <p is defined as zero downward as is 9, and £ is the difference between 9 and </>. 
In what follows we assume that the average of (denoted as £) is zero, and that is 
small. 

The equation of motion for a pendulum driven at any angle, Eq. (jZJ), can be rearranged, 
in the form 

9 = "If sin 6 ~ sin(e ~ 9d) COSUJt = F{ ~ 9) + /( ^' t)- (9) 

We have separated the angular acceleration into a time dependent driving term f(9, t) oscil- 
lating at the driving angular frequency to, and a time independent part F(9) corresponding 
to the effect of gravity. Equation (JHJ) can be rewritten in an expansion in £ for small values 
of£(f): 

+ 1 ~ + + m t) + ^(0, m. (io) 

We keep only the largest rapidly varying terms on each side of Eq. (fTUj) and write 

£~/(0,t) = -i^sm(0-0 d )cosu;fc (11) 
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Note that £ is u 2 larger than so the terms in £ may safely be ignored. Also both on 
the left and F(9) on the right do not oscillate at the driving angular frequency u. 

Equation (JTTJl may be integrated twice, taking the slow motion to be constant on the 
time scale 1/uj: 

£(t) ~ — sin(0 - d ) coswfc (12) 

By averaging Eq. © over the fast component of the motion at angular frequency u, we 
obtain the equation of motion for the slow swinging of the pendulum: 



The rapidly oscillating terms £, £, and f{<p,t) average to zero, and we have 



df 

0^F(0) + ^(0,^. (14) 

We then substitute for F(<ft) and df/dt from Eq. Q and for £(t) from Eq. (jl2j) . Only the 
cos tut terms vary rapidly, and so can be averaged on the longer time scale of the change of 



If we use cos 2 ut = |, we obtain the equation of motion for the angle 

3g . . 9A 2 u 2 



r- 2L Sin< ^ ^-cos(0-# d )sin(0-# d ). (15) 

Equation (fTKJ) describes the slow swinging motion of the driven pendulum. If we use a simple 
trigonometric identity, we obtain another useful form of the equation of motion: 

~ -M. S in - M^! sin [2(0 - 9 d )]. (16) 

The effective torque is the acceleration about one end of the rod multiplied by the 
moment of inertia about the end of the rod: 

r eff (0) = ^0 - ( sin + ^ sin[2(0 - 9 d )]) . (17) 

This torque can be derived from an effective potential energy, V e s{4>), where r eff (0) = 
-dV eS / d<p. 

K ff (0) = ( cos + cos[2(0 - 9 A )\). (18) 

We define the dimensionless critical parameter, R, and the critical angular frequency, u c : 

n 3A 2 u 2 lu 2 , , 

R = —r~r = —2 19a 

AgL u 2 



and rewrite the effective potential from Eq. (JTHj) : 



VM) = |-(cos0 + -cos[2(0-0 d )]). 



(20) 



The first term on the right side of Eq. (j20|) is simply the effect of gravity acting on the center 
of mass of the pendulum. The second term comes from the dynamics of the forced motion, 
and represents the average kinetic energy of the rapidly driven oscillation of the pendulum 
about its center of mass. As the pendulum deviates from the drive angle 9a, the angular 
kinetic energy of the pendulum about its center of mass increases. 

If treated as an effective potential energy, the kinetic energy associated with the driving 
angular frequency stabilizes the slow motion of the inverted pendulum. Figure El shows 
the effective potential as a function of (f> for a driven inverted pendulum with R = 1.75, 
corresponding to the apparatus of Ref. 0. In Fig. 121 we see the gravitational potential 
minimum at <fi = and also the dynamic potential minimum at <fi = 180°. If the drive 
amplitude or angular frequency become too small (R < 1), then the stable equilibrium at 
<p = 180° disappears. (Note that the local maximum near 125 degrees limits the amplitude 
of slow oscillation of this particular case of a driven inverted pendulum.) 

The same physical interpretation of the driven wobble of the pendulum as a stabilizing 
effective torque is seen in Eq. ([17)1 . which includes a stabilizing term proportional to u 2 
whose origin is the kinetic energy of rotation at the driving angular frequency uo. (This term 
is analogous to the centrifugal force in orbital motion which originates from the rotational 
kinetic energy about the center of mass when expressed as a radial equation of motion.) 
Reference U offers additional physical insight into the stability of the inverted pendulum. 

IV. SPECIAL CASES OF THE DRIVE ANGLE 

Let us look at three special cases before moving on to general values of 9&. These three 
examples will provide guidance in interpreting the result for arbitrary values of 9d- 

A. Drive angle zero degrees 

If 9a = 0, we obtain the equation of motion for the slow oscillation <p(t) from Eq. (fT5|) : 




(21) 
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For small angles ~ 0, we can take cos0 ~ 1 and sin0 ~ 0: 

0+^(l + ^)0 = + u, p 2 0-O. (22) 

The pendulum oscillates slowly about = with an angular frequency u p , which is the 
square root of the coefficient of the term in 0(t): 



3A 2 Ld 2 / UJ 2 

u p = u}q\\1+ ^ L =uj Jl + — = uJoVl + R. (23) 



We have used ujq = a/3^/2L, which is the angular frequency of the undriven pendulum. 
Driving the pendulum increases its frequency (decreases its period). 

B. Drive angle 180 degrees 

If 9d = 180° = 7r, we let 5 = — tt and obtain the driven inverted pendulum from Eq. (|15jl: 

5+ (^Z^ cos ^2Z^ sm5 = - (24) 

Small 5 allows us to make the approximations cos 5 ~ 1 and sin 5 ~ 5: 
which gives small oscillations with an angular frequency 



'3A 2 u 2 lu 2 



We obtain stable small oscillations only if R > 1 or equivalently ui > u c . A stable driven 
inverted pendulum has a lower frequency (longer period) of slow motion than a free pendulum 
with the same geometry. 

The traditional treatment of the inverted pendulum starting with Eq. (JJJ) is reviewed in 
the Appendix. The solution of the linearized equation of motion for 9d = 180° gives the 
Mathieu functions. 9 Only for a limited range of drive amplitude and angular frequency do we 
obtain stable oscillations of the driven inverted pendulum. The range of stability is usually 
displayed in terms of Mathieu parameters (see Fig. ITU}). 

Figure H3 displays the regions of stability of the Mathieu functions plotted in terms of 
the pendulum drive parameters A and uj with L = 0.25 m.— The lower smooth curve in 



Fig. El is approximately Auo = 1.81 m/s. Note that we should avoid large driving amplitudes 
(A > 0.08 m) for a 25 cm rod. The effective potential technique of Ref. Q is not valid for 
large excursions from equilibrium, so does not address an upper stability limit. A numerical 
solution of Eq. (JJJ) is needed to explore shorter rods (or large driving amplitudes) together 
with slow oscillation angular amplitudes larger than the linear approximation supports. 

A full numerical solution of Eq. (|7j). including ad hoc frictional damping, for 9d = 180° is 
shown in Fig. |1] Constant friction gives a rapid linear decrease in the oscillation amplitude, 
while viscous damping gives the familiar under-damped oscillation with an exponentially de- 
creasing amplitude. Our demonstration pendulums seem to be better described by constant 
friction. 

Figure shows interesting, non-oscillatory behavior which appears in numerical solutions 
for very short rods. This possibly chaotic non-oscillatory behavior of the driven pendulum 
suggests other studies of overdriven parametric systems which can be simply related to a 
driven mechanical pendulum. Previous studie a 1Ql11 could easily be extended to arbitrary 
drive angles using the analysis in Sec. El 

C. Drive angle 90 degrees 

For a drive angle 6d = 90° = tt/2, the equation of motion for slow oscillations, Eq. (|15jl. 
reduces to a simple form: 



Equation (J28j) gives small oscillations near = only if u < oj c or R < 1. However, if 
uj > uj c (R > 1), the term in parentheses in Eq. (J28|) is negative, and there cannot be stable 
oscillation near = 0. To find the new location of the stable equilibrium angle (0o), we take 
= 5 + 0o i n Eq. (|2"7|) and write 




(27) 



For small angles (cos 




(28) 



5+ -f-[l - Rcos(5 + O )] sin(<5 + O ) = 0. 



(29) 



We obtain the equilibrium angle 0o by taking 5 = and 5 = 0: 



(1 — i?cos0 o ) sin O = 0. 



(30) 
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The solution O = 0, which we have already seen, is stable only if u < u c (R < 1). 
Equation (jSOjl also has a second solution valid for R > 1: 



1- #008 00 = 0. (31) 

The equilibrium angle 0o is nonzero and finite for R > 1: 

0o = cos" 1 i. (32) 

Observe that < O < n/2 for _R > 1 (or w > u c ). 

Near equilibrium this second solution should reduce Eq. (}29|) to an approximate harmonic 
oscillator: 

5 + u 2 p 8 ~ 0. (33) 
Equation ()33|) gives a stable equilibrium only if the coefficient of 5 is positive. 

^ { [1 - R cos(5 + O )] sin(5 + O ) } 5=Q > 0, (34) 

which yields 

[1 — R cos O ] cos 0o + -R sin 2 O > 0. (35) 

The requirement for equilibrium at 0o in Eq. (JBTlj) makes the first term in Eq. (j33j) zero. 
Hence 

J Rsin 2 o >O, (36) 

which holds for all < 0o < 7r/2 (lu > uj c or R > 1). 

Figure|^a) shows the effective potential for R = 1.75 with a stable minimum at 0o = 55°. 
The effective potential for R = 0.75 < 1 is shown in Fig. H2b), where the only minimum in 
the potential is at O = 0. 

The horizontally driven pendulum (9d = 90°) is stable hanging down (0 O = 0) for R < 1, 
and stable near 0o = cos _1 (l/-R) for R > 1. 

V. THE DRIVEN PENDULUM AT ANY ANGLE 

We start with the equation of motion (jlfij) for a general drive angle 9d- 

0+g(sin0 + |sin[2(0-^)]) ~ 0. (37) 
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Equilibrium occurs when = 0o with = 0: 



sin 0o + — sin [2 (0o 



e d )} = 0. 



(38) 



As we have seen, the equilibrium is stable if 



d 

rf0 



(sin0 H sin [2(0 



0*))) 



0=00 



>0, 



(39) 



which gives 



cos 0o + R cos [2 (0o 



e d )} > o. 



(40) 



Small oscillations about the equilibrium angle 0o have the angular frequency u p , where 



The condition for equilibrium in Eq. (|38|) can be easily solved numerically for 9 d as a 
(multivalued) function of 0o- However, we wish to know the stable angle of oscillation O as 
a function of the drive angle 9 d . Simply finding the solutions of Eq. (|3*H|) is not sufficient. 
The zeroes must correspond to real angles and to stable equilibriums as defined by Eq. fl40l) . 
Also, the ideal problem shown in Fig. ^ has a two-fold ambiguity of drive angles at 9 d and 
9(i + 180°. Real lab or demonstration apparatus will limit the motion of the rod so that we 
cannot observe oscillations for \9 d — 0o| > 90°. 

Figure Ufa) shows that a driven pendulum^ with R = 1.75 has a range of drive angles for 
which there are no (nearby) stable equilibrium angles. For all values of R > 2, there is a 
stable equilibrium 0o for all drive angles 9 d - The critical case R = 2 is shown in Fig. Efb). 

For values of R > 2, the curve representing (f)o(9 d ) becomes smooth, tending to a straight 
line as R —>■ oo. (But we must beware of an instability appearing at larger drive parameters, 
equivalent to large R. There are practical limits to how large R can be to give a stable 
driven pendulum in a physical system.) Values of R less than unity gives stable oscillations 
with 9d 7^ 0o over the driving angle range < 9 d < 90° as shown in Fig. |H1 

VI. COMPARING THEORY WITH EXPERIMENT FOR ANY DRIVE ANGLE 

The desire for more degrees of freedom on a limited budget has led us to the driven 
pendulum demonstration using a small hand-held variable speed jig saw (see Fig. EJ). The 
conversion of the jig saw parameters to the units used here gives the corresponding values 




(41) 
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of u and A: < uo < 330rad/s and A = 0.0089 m. The L = 20 cm pendulum gives the 
expected maximum R = 3.3, well into the regime where the driven pendulum is stable at 
all driving angles. 

The value of R can be measured using Eq. (J32j) for R > 1. Hold the driving saw blade 
horizontally and measure the angle to which the driven pendulum rises. At low driving 
speeds, it will oscillate about the downward direction. Once u > u c , the pendulum will 
slowly rise with increasing drive speed until a maximum angle is reached. Any angle above 
60° will correspond to R > 2, which then gives stable driven oscillations at any driving angle. 
The 20 cm pendulum has a measured maximum 0o(9O°) ~ 72°, corresponding to R ~ 3.2, 
in good agreement with our expectation. 

Figure El shows the 20 cm pendulum driven with 6^ = 135° and maximum frequency 
(R ~ 3.2), having being damped into a steady position near 4>o(6d = 135°) ~ 118°. Shorter 
pendulum rods allows the exploration of large R, where chaotic motion of the type shown 
in Fig. 03 may be observed. From simulations and experiments with shorter rods, we find 
that stable driven oscillations do not extend to very large values of R. The unstable system 
shown in Fig. 03 has a value of R of only 14.5. The upper limit on the amplitude from the 
Mathieu function analysis described in the Appendix and shown in Fig. El does not reliably 
predict the onset of instability, because the Mathieu equation only applies for small drive 
amplitudes. The lower limit from the Appendix corresponds to R = 1. Establishing the 
correct upper limit on R for stable oscillation requires further study. 

VII. CONCLUSIONS 

The pendulum driven at any angle from the vertical has been studied using the effective 
potential method of Landau and Lifshitz 1 for rapid driving angular frequency (u ^> uj p ). 
Numerical simulations of the full equation of motion (JI|) generally confirm the results of 
the simplified model for uj ^> uj p . The pendulum will oscillate about the equilibrium angle 
0o defined by Eqs. (J3~%|) and (f4*U|). The angular frequency of small oscillations about 
equilibrium will occur at u p given by Eq. (j4"Tj) . 

The general behavior of the driven pendulum at any driving angle is summarized by the 
parameter R defined in Eq. fll9j) . For R < 1, there no stable inverted oscillations near 
180°. There are stable oscillations with <bn > for all < 9d < 90° and lim <bn = 0. 
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For 1 < R < 2 there are stable inverted oscillations near 180°. Some angles in the range 
90° <6 d < 180° are not stable for \6 d - O | < 90°. Also O (9O°) > 0. For R > 2, there are 
stable oscillations for |0 O — 9 d \ < 90° for all 6 d . 

All parameters of the driven pendulum are accessible over interesting ranges with simple 
and inexpensive apparatus. Even g can be reduced by tipping the plane of oscillation. We 
have explored more parameter sets than reported here and hope that others will explore and 
report their own variations of the driven pendulum at angles other than 0° or 180°. 



APPENDIX A: APPROXIMATE ANALYTIC SOLUTION OF THE INVERTED 
PENDULUM 



Section II V Bl introduced the solution of the inverted pendulum based on the effective 
potential approach of Ref.^. This system also can be easily solved by linearizing the equation 
of motion as in Ref. 

For small angles the equation of motion for the driven inverted pendulum, Eq. (J7J), sim- 
plifies to 

6+(—coaut-^L)5~0. (Al) 
Equation (|A1J) has the form of Mathieu's differential equation 9 : 

d 2 v 

-4 + (a - 2q cos 2z)y = 0. (A2) 
dz 2 

If we make the substitution z = ut/2, Eq. (jA2|) can be expressed in the Mathieu form, 

d 2 5 6o 6A . . . 

The Mathieu parameters can be found by comparing Eqs. (|A2|) and (jA3|) and are q = —3A/L 
and a = —6g/(Lu 2 ). 

The general solutions of Mathieu's differential equation are expressed by the even and 
odd Mathieu functions. For a portion of the parameter space in (a, q), the Mathieu functions 
are real and periodic, corresponding to stable equilibrium of the inverted driven pendulum. 
Outside of this portion of the parameter space, the functions are complex and divergent, 
corresponding to unstable oscillations. 

The region of parameter space that gives stable oscillations is defined through the Mathieu 
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parameters a and q.— 

,2 | ? |3 g 4 U | g |5 

128 ~ 1536 ~ 36864 



q 2 \q\ 3 q 4 ll\q\ , is 



a > a (g) = -— + -i 1- H (A5) 

° w 2. 128 2304 v ; 

For the demonstration pendulum of Ref. |s|, we have q = —3A/L = —0.1524, so q 2 <C 1. 
Then, to a good approximation the lower limit from Eq. (jA5|) simplifies to the leading term. 

a = -u7 2> -^2= — ■ (A6) 

If we solve for the minimum angular frequency for stable inverted driven oscillations, we 
obtain the familiar result of Eq. ()19|) : 



w>We = V^- (A7) 

Including the g 4 correction from Eq. (jA5|) increases ui c by less than 1%. The actual apparatus 
does not warrant this level of precision, so we can safely devise a stable driven inverted 
pendulum using 

Figure is commonly used to display stable regions in the Mathieu parameter space 
of (a,q). It is more useful for our purposes to transform to the drive parameter space of 
(A,u,L). The inversion of the first few terms of Eqs. (|A4|) and (|A5J) for ao(q) and a\{q) 
gives lower and upper limits for the drive amplitude A as a function of the drive angular 
frequency uj and the pendulum rod length L. 



64- J4096- 3584 -gr 



A > A {u,L) = |^ ! — (A8) 



A<A 1 {u,L) = -[ -2 + ^/3 + ^). (A9) 



i(- 2 W 3+ B 

These equations were used to generate the limits of stability in the parameter space (A, to) 
for L = 0.25 m shown in Fig. El 
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FIG. 1: The pendulum driven at angle 9^- 



2 V eff (4>) ImgL d = 180° R = 1.75 




FIG. 2: The potential V e s((p) for a driven inverted pendulum [64 = 180°). 
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50 100 150 200 250 300 350 400 

FIG. 3: The stable region in the drive amplitude A and angular frequency u for a driven inverted 
pendulum of length L = 0.25 m. The parameters of the pendulum of Ref. |8| are shown as a point 
at A = 0.0127m and to = 188s" 1 . 




FIG. 4: Motion for 9d = 180° with constant frictional damping. 
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co= 188rad/s 
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A = 0.0127m 
R= 14.5 
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FIG. 5: Over-driven inverted pendulum of the type modeled in Fig. 01 with L = 0.03 m. No stable 
solution is found in the model or is observed experimentally. 



2V cff W/mgL co 2 = 1.75 <o 2 R=1.75 2^ (<t>) ImgL ar=0.75co 2 R = 0.75 




FIG. 6: (a) Minimum of the effective potential V e s{<p) for R = cj 2 / w c = 1-75 at <po ~ 55°; the 
equilibrium at 0° becomes unstable, (b) For R = 0.75 (u < lo c ), the only minimum in V e s(4>) is at 
<j>o = 0. 
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FIG. 7: All equilibriums (dashed) and stable equilibriums (solid line) are shown in (f>Q versus 0d 
for (a) R = 1.75, and (b) R = 2. 
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FIG. 8: All equilibriums (dashed) and stable equilibriums (solid) are shown in 4>q versus Od for 
R = 0.90. 
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FIG. 9: A driven rod pendulum with R = 3.2 and 9^ = 135° displaying stability at <j)Q ~ 118°. 
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